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Abstract 

At the Accelerator Test Facility (ATF) at KEK intrabeam scattering (IBS) is a strong effect for an electron machine. It 
is an effect that couples all dimensions of the beam, and in April 2000, over a short period of time, all dimensions were 
measured as functions of current. In this report we derive a simple relation for the growth rates of emittances due to IBS. 
We apply the theories of Bjorken-Mtingwa, Piwinski, and a formula due to Raubenheimer to the ATF parameters, and find 
that the results all agree (if in Piwinski's formalism we replace 7^/(3 by H). Finally, we compare theory, including the 
effect of potential well bunch lengthening, with the April 2000 measurements, and find reasonably good agreement in the 
energy spread and horizontal emittance dependence on current. The vertical emittance measurement, however, implies 
that either: there is error in the measurement (equivalent to an introduction of 0.6% x-y coupling error), or the effect of 
intrabeam scattering is stronger than predicted (35% stronger in growth rates). 
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1 INTRODUCTION 



In future e+e- linear colliders, such as the JLC/NLC, 
damping rings are needed to generate beams of intense 
bunches with very low emittances. The Accelerator Test 
Facility (ATF)[|ll at KEK is a prototype for such damping 
rings. In April 2000 the single bunch energy spread, bunch 
length, and horizontal and vertical emittances of the beam 
in the ATP were all measured as functions of current[||],[|||. 
One surprising outcome was that, at the design current, the 
vertical emittance appeared to have grown by a factor of 
3 over the zero-current result. A question with important 
implications for the JLC/NLC is: Is this growth real, or 
is it measurement error? And if real, is it consistent with 
expected physical effects, in particular, with the theory of 
intrabeam scattering (IBS). 

IBS is an important research topic for many present and 
future low-emittance storage rings, and the ATF is an ideal 
machine for studying this topic. In the ATF as it is now, 
running below design energy and with the wigglers turned 
off, IBS is relatively strong for an electron machine. It is 
an effect that couples all dimensions of the beam, and at 
the ATF all beam dimensions can be measured. A unique 
feature of the ATF is that the beam energy spread, an espe- 
cially important parameter in IBS theory, can be measured 
to an accuracy of a few percent. The bunch length measure- 
ment is important since at the ATF potential well bunch 
lengthening is significant[|||. Evidence that we are truly 
seeing IBS at the ATF include (see also Ref. [||): (1) when 
moving onto the coupling resonance, the normally large en- 
ergy spread growth with current becomes negligibly small; 
(2) if we decrease the vertical emittance using dispersion 
correction, the energy spread increases. 

Calculations of IBS tend to use the equations of 
Piwinski[||] (P) or of Bjorken and Mtingwa[||] (B-M). Both 
approaches solve the local, two-particle Coulomb scatter- 
ing problem under certain assumptions, but the results ap- 
pear to be different. The B-M result is thought to be the 
more accurate of the two, with the difference to the P re- 
sult noticeable when applied to very low emittance storage 
rings iQ]. Another, simpler formulation is due to Rauben- 
heimer (R)[^. Also found in the literature is a more com- 
plicated result that allows for x-y coupling[^, and a recent 
formulation that includes effects of the impedance[p^. An 
optics computer program that solves IBS, using the B-M 
equations, is SAD[|ll[|. 

Calculations of IBS tend to be applied to proton or heavy 
ion storage rings, where effects of IBS are normally more 
pronounced. Examples of comparisons of IBS theory with 



measurement can be found for proton[|12||,[|l 3|] and electron 
machines [p4|] , [p3|] . In such reports, although good agree- 
ment is often found, the comparison/agreement is usually 
not complete (e.g. in Ref. [ p^ growth rates agree reason- 
ably well in the longitudinal and horizontal, but completely 
disagree in the vertical) and/or a fitting or "fudge" factor is 
needed to get agreement (e.g. Ref. [Q). 

In the present report we briefly describe IBS calcula- 
tions, and derive a theorem concerning the relative vertical 
to horizontal IBS emittance growths in electron machines. 
We then compare the results of the P, B-M, and R methods, 
when applied to the ATF parameters. Finally, we compare 
the IBS growth in all beam dimensions, including the effect 
of potential well bunch lengthening, for the B-M calcula- 
tion and the ATF data of April 2000. 

Note that this is a revised version of the original re- 
port. After correcting for a \/2 typo found in B-M, and 
after more carefully considering the Coulomb log factor, 
the agreement between measurement and theory has im- 
proved. 

2 IBS CALCULATIONS 

We begin by sketching the general method of calculating 
the effect of IBS in a storage ring (see, e.g. Ref. [^). Let 
us first assume that there is no x-y coupling. 

Let us consider the IBS growth rates in energy p, in the 
horizontal x, and in the vertical y to be defined as 
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Here Cp is the rms (relative) energy spread, Cx the horizon- 
tal emittance, and ey the vertical emittance. In general, the 
growth rates are given in both P and B-M theories in the 
form (for details, see Refs. 
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where subscript i stands for p, x, or y. The functions fi are 
integrals that depend on beam parameters, such as energy 
and phase space density, and lattice properties, including 
dispersion (y dispersion, though not originally in B-M, can 
be added in the same manner as x dispersion); the brackets 
mean that the quantity is averaged over the ring. 
From the l/Ti we obtain the steady-state properties for 
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' We believe that the right hand side of Eq. 4. 1 7 in B-M (with cr^ equal 
to our \/2(Tp) should be divided by \/2, in agreement with the recent 
derivation of Ref. nid]. 
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machines with radiation damping: 
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where subscript represents the beam property due to syn- 
chrotron radiation alone, i.e. in the absence of IBS, and the 
Ti are synchrotron radiation damping times. These are 3 
coupled equations since all 3 IBS rise times depend on e^., 
Ey, and dp. Note that a 4th equation, the relation between 
bunch length as and CTp, is also implied; generally this is 
taken to be the nominal (zero current) relation. 

The best way to solve Eqs. ^ is to convert them into 
3 coupled differential equations, such as is done in e.g. 
Ref. Jl5[], and solve for the asymptotic values. For exam- 
ple, the equation for Cy becomes 



(ky 

dt 



{iy - £yo) 



(4) 



and there are corresponding equations for and Cp. 
Note that: 

• For weak coupling, we add the term —nex, with k the 
coupling factor, into the parenthesis of the Ey differen- 
tial equation, Eq. 0. 

• A conspicuous difference between the P and B-M re- 
sults is their dependence on dispersion ry: for P the 
fi depend on it only through 77^; for B-M, through 
[77' + P'ri/{2/3)] and the dispersion invariant H = 
777^ + 2ariri' + firj' , with a, (3, 7 Twiss parameters. 

• At the ATF, at the highest single bunch currents, there 
is significant potential well bunch lengthening, though 
we are still below the threshold to the microwave 
instability [1]]. We can approximate the bunch length- 
ening effect in our IBS calculations by adding a mul- 
tiplicative factor fpw{I) [I is current], obtained from 
measurements, to the equation relating (Js to ap. 

• The results include a so-called Coulomb log factor, 
of the form \n{bmax/hmin), where bmax, hmm are 
the maximum, minimum impact parameters, quanti- 
ties which are not well defined; typically ln() ~ 20. 
For typical, flat beams we take bmax to be the vertical 
beam size, ay\ bmin = "Tqc^ I {vD = »"o/3x/(7^ex), 
with ro the classical electron radius (2.82 x 10^^'' m), 
c the speed of light, the transverse velocity in the 
rest frame, and 7 the energy factor For the ATF, 
ln() = 16.0. 

• The IBS bunch distributions are not Gaussian, and 
tail particles can be overemphasized in these solu- 
tions. We are interested in core sizes, which we 
estimate by eliminating interactions with collision 
rates less than the synchrotron radiation damping 
rate[p^. We can approximate this in the Coulomb 
log term by letting T^b'^^^^^{\vx\){n) = 1/t, with n 
the particle density in the rest frame [p^; or bmin = 
■^\-Kax<yy(yzl\NcT\ {lixhxf^^, with N the bunch 
population. For the ATF with this cut, ln() — 13.9. 



2. 1 Emittance Growth 

An approximation to Eqs. ^, valid for typical, flat elec- 
tron beams is due to Raubenheimer (R) [H],[[l^: |^ 
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If the vertical emittance is due only to vertical dispersion 
then|] 

eyo ~ Je(Wy)CTpo , (6) 
with J7e the energy damping partition number. We can 
solve Eqs. H^^to obtain the steady-state beam sizes. Note 
that once the vertical orbit — and therefore (Hy) — is set, 
eyo is also determined. 

Following an argument in Ref. [||] we can obtain a rela- 
tion between the expected vertical and horizontal emittance 
growth due to IBS in the presence of random vertical dis- 
persion: The beam momentum in the longitudinal plane is 
much less than in the transverse planes. Therefore, IBS will 
first heat the longitudinal plane; this, in turn, increases the 
transverse emittances through dispersion (through H), like 
synchrotron radiation (SR) does. One difference between 
IBS and SR is that IBS increases the emittance everywhere, 
and SR only in bends. We can write 
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where Jx^y are damping partition numbers, and ();, means 
averaging is only done over the bends. For vertical disper- 
sion due to errors we expect {'Hy)b ~ {T^y)- Therefore, 



)/eyO ^ {T^x)b 



{^x ~ ^xo)l CxO O^x) 
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which, for the ATF is 1.6. If, however, there is only x-y 
coupling, = 1; if there is both vertical dispersion and 
coupling, will be between {Tix) b / {'Hx) and 1. 

2.2 Numerical Comparison 

Let us compare the results of the P, B-M, and R meth- 
ods when applied to the ATF beam parameters and lat- 
tice, with vertical dispersion and no x-y coupling. We 
take: current / = 3.1 mA, energy E — 1.28 GeV, 
Cpo = 5.44 X 10^"*, (Tso — 5.06 mm (for an rf voltage 
of 300 kV), Exo = 1-05 nm, Tp = 20.9 ms, = 18.2 ms, 
and Ty ~ 29.2 ms; fp^ — 1. The ATF circumference 
is 138 m, J, = 1.4, {Px) = 3.9 m, {f3y) = 4.5 m, 
{Vx) = 5.2 cm and {Tix) — 2.9 mm. To generate ver- 
tical dispersion we randomly offset magnets by 15 /im, 
and then calculate the closed orbit using SAD. For our 
seed we find that the rms dispersion {riy)rms = 7.4 mm, 
{Hy) = 17 Tim, and e^o = 6.9 pm (in agreement with 
Eq. For consistency between the methods we here take 
H)=\n[{ay)^^exl{ro {Px))] = 16. 

-Our equation for 1/Tp is twice as large as Eq. 2.3.5 of Ref. |H|. 
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Figure 1: Vertical differential growth rate over 1/2 the 
ATF, forPiwinski (dashes) and Bjorken-Mtingwa (solid). 

Performing the calculations, we find that the growth 
rates in p and x agree well between the two methods; the 
vertical rate, however, does not. Fig. |l] displays the vertical 
differential IBS growth rate 6{1/Ty ), over half the ring (the 
periodicity is 2), as obtained by the two methods (dashes 
for P, solid for B-M). The IBS growth rate 1/Ty is the av- 
erage value of this function. We see that the P curve is 
enveloped by the B-M curve; the average P result is 25% 
less than that of B-M. 




Figure 2: Differential growth rates over 1/2 the ATF, 
as obtained by modified Piwinski (dashes) and Bjorken- 
Mtingwa (solid). 

From the arguments of Sec. 2. 1 , we might expect that we 
can improve the P calculation if we replace t]'^ y/ Px.y in the 
formulation by Hx.y Doing this we find that, indeed, the 
three differential growth rates now agree reasonably well 
with the B-M results (see Fig. As for the averages, the P 
results are all systematically low, by 6%. According to the 
B-M method 1/Tp = 27.0 s-\ l/T^ = 26.0 s^^, 1/Ty = 
19.4 s-i; ap/apo = 1.52, eJeM = 1-90, ey/cyo = 2.30. 
The emittance ratio of Eq. || is = 1.44, close to the 
expected 1.6. 

The dots in Fig. ^,c give the differential rates corre- 
sponding to Eq. H, and we see that the agreement also is 
good. The growth rates in (p,x,y) are (27.0,26.4,19.3) s^^, 
the relative growths in {(Jp,ex,ey) are (1.51,1.92,2.29). 



3 COMPARISON WITH 
MEASUREMENT 

The parameters ap, as, e^, and ey were measured in the 
ATF as functions of current over a short period of time at 
rf voltage Vc = 300 kV. Energy spread was measured on a 
screen at a dispersive region in the extraction line (Fig. ||a); 
bunch length with a sti'eak camera in the ring (Fig.|b). The 
curves in the plots are fits that give the expected zero cur- 
rent result. Emittances were measured on wire monitors in 
the extraction line (the symbols in Fig. ^-c; note that the 
symbols in Fig. ^ reproduce the fits to the data of Fig. ||). 
Unfortunately, we do not have error bars for this data. In x, 
nevertheless, we expect the errors to be small. In y, from 
experience, we expect the random component of errors to 
be 5-10%. As for the systematic component, it is conceiv- 
able that it is not small, since Cy is small and it only takes a 
small amount of roll or dispersion in the extraction line to 
significantly affect the measurement result. 

We see that appears to grow by ^ 85% by / = 3 mA; 
Ey begins at about 1.0-1.2% of exo, and then grows to 
about 3% of Exo, implying that = 1.8-2.4. If we are 
vertical dispersion dominated, with {riy)rms — 10 mm 
and eyo w .012ea;o, then the data nearly satisfies Eq. ^, 

w 1.6. However, normally, after dispersion correction, 
the residual dispersion at the ATF is kept to {rjy)rms — 3- 
5 mm. On the other hand, if we are coupling dominated we 
see that w 1 is not well satisfied by the data. 
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Figure 3: Measurements of energy spread (a) and bunch 
length (b), with K = 300 kV. 

Let us compare B-M calculations with the data. Here we 
take fpw as given by the measurements, and take ln() = 
14. At / = 3 mA we adjust e^o until the calculated ap 
agrees with the measurement. In Fig. ^ we give examples: 
(1) with vertical dispersion only, with {rjy)rms — 7.0 mm 
and e^o = 6.3 pm (solid); (2) coupling dominated with 



3 mm and ej,o — 8.7 pm (dashes); (3) increas- 



ing the strength of IBS by increasing ln() by 35%: i.e. let- 
ting ln() = 19, for the coupling dominated example with 
{Vy)rms = 3 mm and e^o = 14.7 pm (dotdash); (4) same 
as Ex. 2 but assuming a small amount of ey measurement 
error, i.e. adding 0.6% x-y coupling error (the dots). 

We see that, for all examples, ap{I) agrees well with 
measurement, and ex{I) agrees reasonably well also. For 
general agreement for ey (/), we need either a small amount 
of measurement error (e.g. 0.6% x-y coupling measure- 
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Figure 4: ATF measurement data (symbols) and IBS theory 
fits (the curves). The symbols in (a) give the smooth curve 
fits to the measured data of Fig. |[ 



ment error), or for IBS to be 35% stronger than expected. 
A difference in ln() of 5 units implies a factor of 150 in the 
argument. Although there is uncertainty in the Coulomb 
log factor, this difference seems larger than we expect the 
uncertainty to be. Note that the expected error in the IBS 
calculation itself, assuming ln() is correct, is also small: 
l/ln() — 5%[p^. And finally, note that even if we can 
account for the offset by e.g. a 0.6% x-y coupling measure- 
ment error, we see from Fig. Q that the slope of the vertical 
emittance dependence on current is still steeper than pre- 
dicted. 
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4 CONCLUSION 

We have derived a simple relation for relative growth 
rates of emittances due to IBS. We have found that for 
the ATF, IBS calculations following Piwinski (with rf/(3 
replaced by Ti), Bjorken-Mtingwa, and a formula due to 
Raubenheimer all agree well (though one needs to be con- 
sistent in choice of Coulomb log factor). 

Comparing the Bjorken-Mtingwa calculations (includ- 
ing the effect of potential well bunch lengthening) with the 
ATF measurements of April 2000, we have found reason- 
ably good agreement in the energy spread and horizontal 
emittance dependence on current. The vertical emittance 
measurement, however, implies that either: there is error 
in the measurement (equivalent to an introduction of 0.6% 
x-y coupling error), or the effect of intrabeam scattering is 
stronger than predicted (35% stronger in growth rates). In 
addition, the slope of the vertical emittance dependence on 
current is steeper than predicted. 

We thank A. Piwinski for help in understanding IBS and 
K. Oide for explaining IBS calculations in SAD. 
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